plotG(vx,vy,char,size,color)

5773ce0c-4001-4508-9f98-bc8ff5942f2b 25 6 9 decimal

&[DATE] &[TIME] - &[FILENAME]

vx vy char size color plotG n vx length : plot vx 1 el vy 1 el char size color augment : i 2 n range plot plot vx i el vy i el char size color augment stack : for plot 41 line :

Thiele Table

1. Thiele Table ... Thiele(X,Y,Δ) => 'K'

X Y Δ Thiele ε 1018 -^: n X rows : i 1 n range M i 1 el Y i el : for j 2 n range i j n range dM M i j 1 - el M j 1 - j 1 - el - : dM dM abs ε < ε M i j el X i el X j 1 - el - dM / : if : 21 line for for i 1 n range c i 1 el X i el : c i 2 el M i i el : 21 line for Δ 0 ≡ c M if 61 line :

data set XY

data 0.5 - 1.0489 01 0.5 0.9529 1 0.9075 1.5 0.8638 2 0.8217 2.5 0.7812 3 0.7422 3.5 0.7048 4 0.6689 4.5 0.6343 5 0.6012 5.5 0.5694 6 0.5389 6.5 0.5096 7 0.4816 7.5 0.4548 8 0.4292 8.5 0.4047 9 0.3813 9.5 0.3589 10 0.3376 10.5 0.3172 11 0.2978 11.5 0.2755 12 0.2546 12.5 0.2354 13 0.2176 13.5 0.2012 14 0.186 14.5 0.172 15 0.159 15.5 0.147 16 0.1359 16.5 0.1256 17 0.1162 17.5 0.1074 18 0.0993 18.5 0.0918 19 0.0849 19.5 0.0785 20 0.0726 25 0.032 30 0.015 32 0.011 35 0.006705 40 0.003144 45 0.001536 47 0.001165 50 0.000798 52 0.00062 55 0.000432 60 0.000233 61 0.000205 71 0.00005242 85 0.00000568 562 mat :

ListCollect(data,select)

data select ListCollect collects index list ct 1 : i select a ct el data i 1 el : b ct el data i 2 el : ct ct 1 + : 31 line for a b augment 41 line :

Cfr(j,K,X,x)

j K X x Cfr Optimiz Symbolic j rows(K)-2 N j 1 + el K j 1 + el x X j 1 + el - K j 2 + el / + : c j el N j 1 + el : a j el K j el x X j el - c j el / + : jj j 1 - (1 range c jj el a jj 1 + el : a jj el K jj el x X jj el - c jj el / + : 21 line for a 1 el 71 line :

Error Stem Plot

Error stem plot for QuickPlot unicolor data Error stemX 1 data rows range : i 1 stemX rows range stem i el data i 1 el data i 2 el data i 1 el 022 mat : for i 1 stem rows 1 - range stem i el stem i el stem i 1 + (1 el 1 el 012 mat stack : for V stem 1 el : j 2 stem rows 1 - range V V stem j el stack : for V V data data rows row stack : 71 line :

Select the support points select 2691215192317 mat :

XY data select ListCollect :

Plot data & support points

X XY 1 col : Y XY 2 col : Kvec X Y 0 Thiele : 31 line

Pts i 1 data rows range p i el . : for data 1 col data 2 col p augment 21 line :

X Kvec 1 col : K Kvec 2 col : 21 line

support X Y + 15 blue plotG :

Support points XY 012 0.8217 3.5 0.7048 5 0.6012 6.5 0.5096 8.5 0.4047 10.5 0.3172 72 mat

Thiele 'K' convergents Kvec 012 11.21705 - 3.5 2.346248 - 5 20.137099 6.5 0.8999 8.5 3.136171 10.5 2.051205 72 mat

Activate Thiele

j Kvec rows 2 - :

j 5

x g j K X x Cfr :

Relative residuals for the "ErrorPlot" Residuals sub data 12412 submatrix : i 1 sub rows range δ i el 1 sub i 2 el sub i 1 el g / - eval : for sub 1 col δ augment 31 line :

Pts x g support 3 1 sys

Relative residuals explorer

Residuals Error Residuals 2 1 sys

Export an easy 'J(x)' fraction
evaluate coefficients from symbolic ... 6 decimals adequate x J 1.656485 180.575067 x 43.69696 + 1299.569741 x 5.616953 + 0.000131 x 5.116572 - / + / + / - :

1.604858 174.727169 x 43.473225 + 1286.481982 x 5.24242 + 0.000159 x 5.14073 - (/ + (/ + (/ -

Optimiz  ... NONE x g confrac x convert maple 249304002720682613029772340587685402131155343362927792238765828883903950000000 / 337315137813063961200144298337459873443111874881367749361144090183307176827719305248324572620860463693919042711192695819601415781914260482000000000000 / x 2085897463274697688895878614842201111959009173348707853502567124096407801900959654435519335754907436791374774798119919359305013110642132904636926782891001375131557050165632299580308653759683173108765519383469960000 / + 19176993664769020159021959545938711284739327100713607201458313342258977639520650794392998169946505898949197469798188972490062628669599390029420951490653886801414864882647627404376934555600574364646605799920762268278179849786835525995668700546438031016065567686877765834669464061472292361 / x 60591785807636540652835674635269653629581891212275001552712942536329189652240050081324356057220182306444911185645750024626913436294875362135837976625523522352940835582879675144809621573471286900127248816901155797923684286325933492322527289757050918001256352543068383967162378579566250795389192023190144087454651768735727311354361693300596837837754001749694508139350285431950834918598790123259105448466966847543 / + 142355177841333885215163420703151546413011681818372535618136287850673461465628238851219935325079994690691703188803789405143362312475937326953100517861291945751821779068697973609668945097529474658645921508400954292460851731741064973996787218657704000000000000000000896162987411759927818203984963715939540606801528170920472096458697601176054786636972055723736430037319288663661525223955141966161689360278040075996712850883199829948485875311144030325411739570874945635246958883794179466924885388886244558473539339445537775042669449609 / x 30778558351128879859951234135280636394744926103294090493881611607547681849409534514295136113048867835501825401793512763883064738565248159871962968045733779028932207669315538007994679485969697676710637242202942126733888331012933625765770963493721096123045931295322275706 / - (/ + (/ + (/ -

19176993664769020159021959545938711284739327100713607201458313342258977639520650794392998169946505898949197469798188972490062628669599390029420951490653886801414864882647627404376934555600574364646605799920762268278179849786835525995668700546438031016065567686877765834669464061472292361 / 1286.481982

447460047466545858074218793653703001072698796742311501762396101001878350210613946127576255467004887041791434600932284847435956851260431410095343273537771575374530335843420153736757541931136419033509122126861216425826571173062447235231763297587801272443684785021214934798454909161 / 1303.508713